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We consider the problem of quasiparticle poisoning in a nanowire-based realization of a Majorana 
qubit, where a spin-orbit-coupled semiconducting wire is placed on top of a (bulk) superconductor. 
By making use of recent experimental data exhibiting evidence of a low-temperature residual non- 
equilibrium quasiparticle population in superconductors, we show by means of analytical and nu- 
merical calculations that the dephasing time due to the tunneling of quasiparticles into the nanowire 
may be problematically short to allow for qubit manipulation. 



Devising a physical system where to experimentally ob- 
serve for the first time the presence of Majorana fermions 
has become in the last years a serious and stirring chal- 
lenge in the condensed-matter physics community. Apart 
from the importance per se of observing the signature 
a Majorana fermion, the strong hope is to realize a 
Majorana-based qubit, which would offer an intrinsically 
improved protection against decoherence due to the pe- 
culiar delocalized structure of the Majorana state itself. 
For example, in the setup we consider ,1^11^ a topolog- 
ical superconducting state (TSC) in a semiconducting 
nanowire is created, with a mid-gap mode dl^^ at en- 
ergy Eond — 0, whose wavefunction is strongly localized 
at the two ends of the nanowire. Such topological state 
can be induced through the combined effect of s— wave 
pairing, spin-orbit coupling and magnetic field^^. The 
superconducting pairing is inherited most typically via 
proximity effect from a bulk superconductor placed be- 
low the wire (around the wire in some proposals). Even 
without restricting ourselves to this specific setup, su- 
perconductivity is a key ingredient needed in essentially 
all the proposals to produce observable Majorana excita- 
tions in condensed-matter systems. The considerations 
we make in this paper are thus qualitatively valid and rel- 
evant for a wide range of configurations, while the quan- 
titative results are specific to the proximized nanowire 
setup. 



The zero-energy many-body excitation dcnd in the 
TSC can be exploited to store information in an ideally 
dephasing- free qubit. Defining the |0) state as the many- 
body state where the dend is empty, and correspondingly 
|1) = dondl'^)' the subspace spanned by |0) and |f) is a 
degenerate ground-state subspace, which offers intrinsic 
protection against dephasing. However, coherent states 
of the type (|0) -I- e'^^jf )) cannot be prepared, because 
there is no physical coupling which could create such su- 
perposition. Strictly speaking, then, a system where the 
states |0) and |1) differ by fermion parity (occupation of 
a single BCS-hke mode) cannot be used as a qubit. To 
obtain a proper quantum bit, one needs at least two of 
these zero-energy states, that is, four Majorana fermions. 



In such case, there are four degenerate states: 

|00) = |0)end,l ® |0)cnd,2 , 

|10) ^4„d,il00), 
|oi) ^ <„d,2l00) , 
111) ^rfUi4„d,2l00). 

The states |00) and |1I) (possible choice for the qubit) 
share now the same fermion parity, and if we choose them 
as qubit computational states, coherent superpositions 
are possible thanks to the superconducting pairing which 
induces fluctuations in the number of electrons, in jumps 
of two, due to the hopping in and out of Cooper pairs. 

If the superconductor can only exchange Cooper pairs 
and not single, unpaired electrons with the wire, then the 
fermion parity (i.e. the number of electrons modulo two) 
is conserved. This is at the base of the protection these 
systems benefit!^. If instead single electrons could enter 
the nanowire, then the system would be driven out of the 
topological subspace {|00), |11)}, populating |10), |01) or 
some higher-energy states. In the case of a single TSC 
segment, with only one zero-energy mode, the presence of 
unpaired electrons would instead cause errors, causing 
transitions |0) |1) and viceversa. 

There have already been a couple of work&^iSJ point- 
ing out that Majorana-based qubits are prone to stan- 
dard decoherence mechanisms when one allows for single- 
electron tunneling from a generic external (noisy) envi- 
ronment. The specific phenomenon of the possible dis- 
turbing presence of unpaired electrons in the supercon- 
ductor, dubbed "quasiparticle poisoning" , constitutes a 
similar but more subtle problem, since the tunneling elec- 
trons are coming from the superconductor itself, which 
is an indispensable ingredient in these type of propos- 
als. The issue has already been taken into consider- 
ation in the Majorana- fermion communitj'^Hn] but, to 
the best of our knowledge, no explicit estimation has 
been made for a specific setup. Rather, some estimates 
for quasiparticle tunneling rate from the superconductor 
into the wire have been borrowed from some recent ex- 
periments dealing with equilibrium superconductivitjff^. 
In such experiments, like in several other investigations 
beforel^HiSl^ low-temperature quasiparticle poisoning has 
been observed, and evidence has been provided that the 
quasiparticle density does not drop to zero at the small- 
est temperatures, as predicted by theory, but rather sat- 
urates to a constant value. A clear and widely accepted 
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explanation for this phenomenon is still lacking, but the 
observation of excess quasiparticles has been repeatedly 
confirmed. For instance, a very recent experimental in- 
vestigation on quasiparticle kinetics inducing energy re- 
laxation in a transmon qubit has measured poisoning 
times below fi^. 

In an ideal superconductor at zero temperature all the 
electrons are forming Cooper pairs and out of the conden- 
sate no quasiparticles exisll^. At finite but small tem- 
perature, the average density of quasiparticles in a super- 
conductor at equilibrium follows the activation behavior 
{kB = 1) 

nil = 2:.sV27rrAsexp(-As/r) , (1) 

valid at T < Ag, with the normal-state, single-spin 
density of states at the Fermi level, T the temperature 
and Ag the energy gap of the superconductor. The corre- 
sponding average quasiparticle recombination time (life- 
time) is evaluated aJ^^ 

r. = ^f^)'yf e-^Z-oc^, (2) 
Vtt v2As/ V T riqp 

where Tc is the critical temperature of the superconduc- 
tor and To a material dependent, characteristic electron- 
phonon interaction time. Eqs. ([T]) and ^ predict a very 
low quasiparticle density and correspondingly a very long 
quasiparticle lifetime at temperatures T <C Tc. As an- 
ticipated above, experimental data show agreement to 
the exponential behavior Eq. ([2| only for not-too-low 
temperatures. For instance, the work Ref. [12] obtains 
a zero-temperature saturation Tr ~ 2 ms in the quasi- 
particle lifetime for temperatures below 160 mK, and 
a corresponding saturating quasiparticle density n^p ~ 
25 — 55 /Ltm~'^. These large quaisparticles lifetimes of the 
order of ms have been sometimes used in the literature as 
an estimate for the tunneling rate of quasiparticles into 
the TSC nanowire, i.e. an estimate for the qubit lifetime, 
leading to the conclusion that quasiparticle poisoning is 
not a serious issue. 

Such experiments-^, however, analyze the quasiparticle 
lifetime in an isolated superconducting system (two- 
dimensional in the case of Ref. [T^] ) , and not in a hybrid 
structure where two subsystems are put into contact, and 
a subsystem can be poisoned by the other one. 

Here we would like to stress that the rate of tunneling 
into the qubit subsystem (the nanowire), and thus the av- 
erage dephasing time of our Majorana qubit, is not given 
by the quasiparticles lifetime as measured in an isolated 
superconductor. Indeed, for example in the case of our 
SC/TSC junction, quasiparticles in the superconducting 
reservoir can well tunnel into the TSC wire and tunnel 
out again, many times before a recombination event may 
take place after the typical time Tr. Since even a single de- 
tour of the quasiparticle into the qubit part of the system 
could destroy the coherence of the qubit itself, it is impor- 
tant to calculate or measure explicitly the tunneling rate 
of quasiparticles Fqp into the nanowire. To this end, we 



adopt a well-established formalism, already employed in 
earlier calculations for Fqp which were performed for the 
case of Josephson qubits, adapting them to the present 
case of a SC interfaced to a TSC nanowire. We demon- 
strate that, depending on the parameters, the tunneling 
rate of quasiparticles can vary in the range 0.1 — 100 MHz, 
imposing therefore a much more serious constraint on the 
feasibility of error-free qubit manipulation. This becomes 
particularly clear in the last section of this paper, where 
we produce quantitative estimates for the poisoning rate 
in typical configurations, and we compare them with the 
time scales required for adiabatic computations. 

I. CALCULATIONS 

The system we consider is constituted by a bulk 
s— wave superconductor in tunnel-contact with a TSC 
nanowire. Gap magnitudes As and At in the SC and 
in the TSC need not be the same. Rather, the topolog- 
ical gap At is always smaller than Ag, for two reasons: 
(i) the proximity-induced gap is in general smaller than 
the parental gap Ag, depending on the transparency of 
the interface SC/TSC?W; (ij) the topological p— wave gap 
At is only a fraction of the induced s— wave amplitude, 
and strongly depends on the values of magnetic field, 
spin-orbit coupling^, and even on the electron-electron 
interactiorPi^. 

The bulk superconductor is described by a standard 
BCS Hamiltonian i/gc, whose eigen-excitations are 
bogolons /3ka with energy Es{k) = i/^fTAf = 

i^sc-E^s(fc)/3L/3fe- (3) 

k(T 

As mentioned above, the nanowire can be driven into a 
topological superconducting phase by means of the com- 
bined effect of spin-orbit coupling, Zeeman interaction 
and superconducting pairing^l. The topological phase is 
reached by tuning the chemical potential within the gap 
opened by the Zeeman interaction between the two chi- 
rality bands, and singling out in this way a single pair 
of Fermi points. In terms of this single effective de- 
gree of freedom, the original s— wave superconducting 
interaction becomes of p— type. The low-energy effec- 
tive Hamilton ian f or such spinless p— wave superconduc- 
tor then read^^J^ 

H^w = Y ^Skdldk + sgn(fc) (^AT4^tfe + ^Td-kd-k^ , 

k 

. . 

where the d^'s describe the lower-band electrons originat- 
ing from the combined effect of spin-orbit and Zeeman in- 
teraction, with dispersion Sk- After diagonalization, the 
low-energy Hamiltonian of the nanowire is also written 
in terms of bogolons rjk , 

H^w Htsc = Yl ^T^i'^)^kVk , (5) 

k 
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with dispersion E'T{k) = \/s\ + A^. Finally, the two 
subsystems are coupled by tunneling, described by the 
Hamiltonian iJxj 



(6) 
(7) 



The operators Op^o- annihilate an electron in the state 
\p,cr) in the SC reservoir, while the Ck^a-^s are bare- 
electron operators in the nanowire. Switching to the di- 
agonal basis of Zeeman and Rashba in the nanowire leads 
to the final expression ([?]) written in terms of the effec- 
tive spinless lower-band electron operators dk introduced 
above. Here the tunneling amplitudes tkpa differ from 

the bare-electron tunneling amplitudes i^^*, since they 
describe the hopping into the effective spinless modes 
dk , and incorporate the spin-dependent factors which de- 
scribe the mixing of degrees of freedom due to spin-orbit 
and Zeeman interaction. 



II. ESTIMATION WITHOUT 
ENVIRONMENTAL P{E) THEORY 

We first start with the case where the only relevant 
degrees of freedom are those related to the electronic 
quasiparticle tunneling through the SC/TSC junction. 
In a more refined theory the event of a quasiparticle tun- 
neling through the junction is influenced by the charge 
dynamics in the environment around the junction itself. 
This approach, the so-called ^^environmental P{E) the- 
ory"^, will be considered separately below. 

In order to estimate the rate Fqp of tunnel events from 
the superconductor to the nanowire, we start with a 
Fermi's golden rule approach, along the lines of Refs. 
I27j . and evaluate Fqp by averaging the matrix elements 
of the tunnel Hamiltonian over initial and final configu- 
rations with the appropriate thermal occupation factors: 

- Y E K/I^tN) \'S (Ef - E,) f{E,)[l - f{Ef)] . 

(8) 

The initial and final states \i) and |/) are eigenstates 
of Hq — Hsc + Htsg- The TSC state in the nanowire 
is induced by proximity effect, microscopically described 
by the same Hamiltonian that we are considering 
now. There is, however, no inconsistency, since the first 
order contribution Eq. ([s]) does not take into account 
the Cooper-pair hopping, which is assumed to be already 
implicitly included in i/xsc- 

We are interested in calculating matrix elements of the 
type 

{f\HT\t) ^ {nk^l,np^O\HT\nk^O,np^l) , (9) 



where we have indicated by jrifc = 0, rip = 1) the prod- 
uct state with the TSC in its ground state (zero quasi- 
particles) and with one excess quasiparticle in the state 
\p) in the bulk superconductor. Correspondingly, \nk — 
1, Tip = 0) describes the state where the bulk SC is in its 
ground state, and one quasiparticle \k) is now present in 
the nanowire. The matrix elements of the above equation 
can be evaluated by using the Bogoliubov transformation 
which diagonalizes the BCS Hamiltonian, 



'^pi^p.a ^~ ^'^pi^—p,<T 
Cipher '^pl^p^cr ^'^pl^—p^<7 



(10) 



Here u and v are the usual particle-like and hole-like co- 
herence factors, and a = —a. An analogous transforma- 
tion can be applied to the dj. operators in the nanowire, 
with corresponding u/v coefficient^^. The explicit ex- 
pression for the coherence factors is (we now denote them 
by us,T and 113, t in order to make clear to which subsys- 
tem they refer to) 



ul{E),vl{E) 



1 ± 



A2 



E 



(11) 



and similarly for ut{E) and vt:{E). After this step the 
Hamiltonian formally describes tunneling of quasi- 
particles /3p.o- and rjk across the junction. The insertion 
of the Bogoliubov transformation into the Fermi's golden 
rule produces the formula 

Fqp=^ H \tkpa\^[u{Ep)u{Ek) - v{Ep)viEk)f 



k,p,(7 



2Tl\t\ 



xr^{Ep)[l-f{Ek)] 6{Ek-Ep) 
Y,C{Ek, Ep)r^{Ep) 6 {Ek - Ep) 

k,p 



(12) 



Note that we have added a superscript to the Fermi oc- 
cupancy factor in the bulk superconductor, to emphasize 
that its quasiparticles follow a non-equilibrium distribu- 
tion, corresponding to the observed excess quasiparticle 
density. Nevertheless, /"'^'^(E') is still assumed to ex- 
hibit a sharp jump at -E = Ag. Further considerations 
about /"°'' are developed in the following sections. As- 
suming a weak energy- and momentum-dependence of 
the tunneling amplitude for energies close to the Fermi 
energy, we have extracted ifc.p out of the summation and 
replaced it with an average squared tunneling amplitude 
l^p. Further, we have made use of the low-temperature 
assumption to discard the term P^°'^{Ep)f{Ek), since, as 
we will recall later in the paper. Ax is typically only a 
fraction of As, and T <C (As — At). Finally, the func- 

iionC{E,E') EE [us{E)ut{E')-vs{E)v^{E')Y has been 
introduced for brevity. 

Converting the sum into integral and using the delta- 
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function constrain gives us 




AsAi 



^{E^-Ai){E^-A 



r^{E) . (13) 



The superconducting density of states J^(-E) in the two 
subsystems is related to the normal-state density of states 
v'^{E) through the expression (j = S, T) 



KiE) 



E 



E^ 



(14) 



The dimensionless tunneling conductance is defined 
as /i/(e^i?T) = Rq/Rt, with Rq the quantum of resis- 
tance and i?T the normal-state resistance of the tunnel 
junction, determined by the formula 



h 



4^6^^ I- 12 



k.p 



47re2|t|2i.^^(0)i/^(0) 



(15) 

We would like to connect at this point the final expression 
for Fqp to the non-equilibrium density of quasiparticles, 
and use the experimentally measured values as an input 
for the calculation. Using the assumption that /"'''^(i?) 
is exponentially peaked at E' = Ag, we can approximate 



all the well-behaved factors in the integrand of Eq. ( 13 ) 
by their value at i? = As- Then, recalling the connection 
between the non-equilibrium quasiparticle density in the 
superconductor and the non-equilibrium Fermi distribu- 
tion (from now on we will simply write n^'P = n^^^), 



noq 



2 / dE iysiE)r^iE) 



(16) 



we can extract a factor n^'^'^ from the integral Eq. (13 1, 
and relate Fap dire ctly to the observed excess quasipar- 
ticle densit> 'i5 | 24 | 25 t 



3T 



,ncq 



As 



As 



(17) 



The square root factor is of order unity for typical values 
of At. Plugging in at this point the experimental values 
for n"°'^ ^ 10//im^ (Ref. [H]), normal-state density of 
states z/g ^ 10^/(^m^ • K), and typical values for gT in 
phase-qubit experiments (i?T ~ 10^ O o f/T ~ 10^ in 
Refs. |28)29j ). we obtain an estimation for Fqp of the 
order of - 10 MHz. 

If instead one has higher tunnel resistances (gx ^ 1 
in Ref. [SU] and gr ~ 10 in Refs. |31I320 . then the rate 
can be largely suppressed. Simply increasing the tunnel 
resistance however does not constitute a valid strategy in 
our situation, because also the tunneling of Cooper pairs 
would be reduced in that case, lowering the topological 
gap in the nanowirc. For a more detailed discussion of 
this point, we refer the reader to the final section of the 
paper. 



III. ESTIMATION WITH ENVIRONMENTAL 
P{E) THEORY 

We now take into account the fact that the tunneling 
probability for a quasiparticle is influenced by the cou- 
pling with the surrounding environment, by making use 
of the environmental P{E) theory'^. This amounts to 
starting with the modified tunneling Hamiltonian 



-tip 



H.c. 



(18) 



where the charge displacement operators e^"^"^ act on the 
electrical circuit degrees of freedom (environment), and 
describe the transfer of a ±e charge through the SC/TSC 
junction in a tunneling event. Here, if is the conjugate 
coordinate to the charge q, with commutation relation 
[q, if] — ie, and gives a charge displacement operator ac- 
cording to the relation e'^^'^qe~^'^ = q — e. Rewriting 
Ht in term of Bogoliubov operators, we obtain several 
terms, among which the ones describing the transfer of a 
quasiparticle have the forrrPSl 

(usuxe-"^ - vsVTe'"^^ rjlPp^^ . (19) 
The evaluation of the modified tunnel rate 

J2\{f\HT\i)\'S (Ef - E,) r^iE,)[l f{Ef)] (20) 

now involves also averages over environment degrees of 
freedom, and it requires the calculation of the equilibrium 
correlation function 



,([v(t)-¥.(0)]v(0)> 



(21) 



which in the case of Bogoliubov-quasiparticle tunneling 
must be properly corrected, as explained by Martinis et 
alW^^ and becomes 



Jit) 



X c 



{u^ + t;2)e<^(*)^(o)) - 2™e-<^(*)^("» 



(22) 



Here for sake of brevity we wrote u = ugux and v — fs^x- 
The fluctuation-dissipation theorem provides us with a 
relation between the correlation function J(t) and the 
dissipation in the environment, indirectly described by 
its impedanc^^E 



J(t) = ([^(i)-^(0)](p(0)) 



= 2 



1) 



(23) 



where Zt^uj) is the total environmental impedance, Rq — 
hje^ is the quantum of resistance for single-electron 
charge transfer, and we have assumed T = (while it 
is still necessary to use a finite value of T in the tunnel 
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rate calculations). This description in terms of circuitry 
elements, where the tunnel junction is characterize by 
its capacitance C and tunnel resistance Rt, and the en- 
vironment properties are encoded in its impedance, is 
summarized in Fig. [l] Due to the presence of the delta 
function in the summations of Eqs. (Isl) and (20 1, what 



we finally need is the Fourier transform of J{t) (and J{t) 
respectively), usually named P{E): 



dt 

2TTh 



.J{t)\Et/h 



e"-'e~-'^ = P{E) . 



(24) 



In terms of such function, the tunneling rate in the case 
of electron-environment coupling is expressed as 



dE / dE'y^ 

As •.'At 



(25) 



ys{E)ME')r''{E)[l - I{E')]P{E - E') 



where now P{E^E') may be interpreted as the probabil- 
ity of a tunnel event which involves an energy exchange 
{E — E') between quasiparticle and environmental de- 
grees of freedom (to be precise, the energy E — E' is 
the energy transferred from the tunneling particle to the 
environment) . 

Going back now to the modified correlation function 
J{t), few comments are in order. The first term of the 



rhs of Eq. (22 1 equals (w^ +v'^)e'^^*\ and the same steps 
described above lead to contribution (m^ -I- v^)P{E) = 
C+{E)P{E). The correlator {ip{t)ip(fi)) is not well de- 
fined, due to an infrared divergence, for impedances 
whose real part does not vanish at = 0, see Eq. (23 1. 



This does not constitute a problem for the physical quan- 
tity J(i), since there the diverging static correlation Lp^ is 
subtracted off. The same does not happen with the sec- 
ond term, which instead involves the factor e~^''''^*^'^^°^\ 
not compensated by e~^f{o)v{'^)) ^ Since {ip{t)iplfi)) is pos- 
itively diverging, however, this second term in Eq. ( |22[ ) 
vanishes. In Ref. [S^ this issue is not present, since they 
consider a model where 3fieZt(0) — 0, and the divergence 
is absent. 



A. Single-mode environment 

One can get the simplest model for the environment by 
studying the coupling of the tunnel junction to one single 
environmental mode, which could come from a resonance 
in the lead impedance Z(uj) of from bound states in the 
barrier. Such coupling can be implemented by putting 
an inductor with inductance L into the external circuit. 
Seen from the junction, the impedance Z{u)) = iujL is in 
parallel with the capacitance C of the junction itself, and 
the total impedance read^ 



1 



1 



(w - ie)2] 



(26) 




FIG. 1: Schematic circuit representation of the 

Superconductor/Topological-nanowire system. The interface 
which separates the two subystems acts as a tunnel junction, 
with tunnel resistance -Rt and capacitance C . The inter- 
nal impedance Zs and Znw of the superconductor and the 
nanowire are combined in the text in a single global environ- 
mental impedance Z{uj). An external voltage bias between 
the two sides of the junction can be present. 



with environmental resonance frequency = 1/ ^ EC . 
The infinitesimal imaginary part e is necessary in or- 
der to obtain the correct result for the real p art o f the 
impedance. By taking the limit e — > one ge 



3fie [Zt(c^)] 



2C 



[(5(w -I- Wr ) -I- J(clI — Wr )] 



(27) 



This expression is essentially saying that the resonator 
can both absorb or emit photons (mode quanta) at fre- 
quency wr. In our case, if we identify the environmental 
mode with the only available low-energy excitation in 
the nanowire-superconductor system, i.e. the Majorana 
mode, we obtain a resonance energy fiwR ~ (or energy 
much smaller than all other energy scales) . This situation 
is sketched in Fig. [2] 




FIG. 2: Equivalent circuit for the SC/TSC system in 
the case of a single-mode environment. The environment 
is modeled with a single inductance L with impedance 
Z(lS) = iojL, corresponding to a total impedance 5Re [Zt(a;)] ~ 
+ cjr) + — i^r)] (plotted in the inset) . 

Before inserting this form of the total impedance in the 
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formula Eq. (24) for P{E), let us define the parametei'^SI 

1 



^ Ec_ 

cjR RqC hoj-n. 



(28) 



which compares the single-electron charging energy with 
the environmental mode excitation energy. This param- 
eter determines the size of charge fluctuation^^!! 



coth ( ^ 

4C V 2r 



(29) 



Using the definition of C, the expression for P{E) in the 
single-mode limit becomes 



P{E)^ 



dt 



exp [C (e 



1)] . (30) 



In typical superconducting charge-qubit experiments one 
haa^Si that Eq <C ^wr, that is C ^ 1, and then the ex- 
ternal exponential in Eq. ( 30 ) can be expanded around 
zero. In our case instead, since ~ 0, such simplifica- 
tion is not possible, and we rather expand the internal 
exponential for t <^ h/uj^. The result is: 



Jit) ~ C {-i^^nt - WRt^) 



E, 



{it -f uj-Rt^ 



(31) 



P{E) 



[E-EcY 
iEchu>-R 



exp 

\/ Ait Echius, 



-^5{E-Ec) . (32) 



That is, in first approximation the energy exchange be- 
tween quasiparticles and environment occurs with unit 
probability and is peaked at the charging energy Eq = 
e^/(2C). In the opposite limit C 1 vahd for typical su- 
perconducting qubit^i^l^, one would get instead that the 
energy exchange is peaked at the resonator energy wr, 
and that the probability ^ C for such exchange is very 
small (the most probable event being the tunneling of a 
quasiparticle without energy flow to the environment). 

Plugging now into Eq. (25) the form of P{E) just ob- 
tained, we get 



for Fqp in the case of typical values for C \ pF) 
and Eq ~ 0.1 /ieV <C As, At. Thus, even in this case 
the typical tunneling rate turns out to be Fqp ^ 100 
kHz— 10 MHz, depending on the transparency of the 
tunnel barrier. 



B. Ohmic environment 

We consider now the more realistic case of an Ohmic 
environment, with external impedance Z{u)) — R and 
total impedance 



1 

iLoC + l/R 



1 



1 



9 [1 + {lu/lucY 



where we have introduced the dimensionless environmen- 
tal conductance g = Rq/R and the frequency 



UJQ 



1 

RXJ 



g_Ec 

TT h 



(35) 



which represents an effective cutoff for the total 
impedance, due to the junction capacitance: at ener- 
gies small compared to huJc the real part of the total 
impedance is essentially given by i?, while for higher 
energies Ztiu) decreases. Such behavior is shown in 
Fig. |3] The P{E) corresponding to this case cannot be 




/iF„ 



[E [E + Ec) - AsAt] 



gr I dE 

/As J[{E + Ecf-A^[E^~A 



(33) 



This is essentially identical to the previous result 
Eq. ( [l3| , with the simple substitution Ep — ^ {Ep + Eq), 
and leads to the low-temperature result 



/iF 



,ncq 



qp 



5t 



An 



Ec 



V As + At + -Ec 



(34) 



Compared to equation (17), the presence of Eq produces 



a negligible modiflcation to the quantitative estimate 



FIG. 3: Same as in Fig. ([2|, but with a different environmen- 
tal impedance. Here the environment is modeled hy Z — R, 
and the resulting Lorentzian total impedance is shown in the 
inset, where uuc ~ 1/{RC) — gEc/irh. 

calculated analytically, but the low-energy and the high- 
energy asymptotic behaviors can be obtained a^^S 



P{E) = { 



e-27/s 1 
r(2/g) E 



TT E 

3^ 



2/s 



7r2 E^ 



for E <t:Ec 



for E^ Ec 



(36) 



Here 7 is the Euler constant and F the gamma-function. 
The behavior of P{E) for intermediate energies has to 
be evaluated numerically. Since we are mostly interested 
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in energy exchanges between the superconductor and the 
topological nanowire of the order of SA = (As — At) ~ 
©(Ag) ^ meV, and since most typically Eq <^ As, we 
are not justified to use the small-energy expansion of 
P{E), and we must rather determine P{E) numerically. 
By taking the derivative of Eq. (l23| and performing a 
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g 


= 1 




= 2 




= 20 




FIG. 4: Behavior of the probability function P{E) as a func- 
tion of the energy exchange measured in units of Ec- For 
large environment resistances (small g) the junction releases 
a typical energy amount of the order of the charging energy. 
For small resistances (large g) the energy which is exchanged 
shrinks to zero, and one recovers a situation with independent 
quasiparticles and junction degrees of freedom. The curves 
have been obtained through numerical integration. 

Fourier transform, one finds that P{E) obeys to the in- 
tegral equation 



EP{E) 



2 



AE' Ike 



E~E' 



P{E') , (37) 



which can be solved for example by iteration. A collec- 
tion of solutions for different values of the parameter g 
is shown in Fig. |4j Qualitatively different behaviors are 
observed in the highly resistive and low-resistive limits. 
Inserting the obtained solution P{E) into Equation ( 25 ) 
we can get the desired estimation for Fqp in this case. 
However, due to the finite energy exchange allowed by 
P{E), the singularities of the two density of states dis- 
tributions can overlap in the integral, and caution must 
be exercised. In particular, one cannot always make the 
simplification adopted to attain Eq. (17), which allowed 
us to single out a factor 

In the case Ec <C 5A the same approximation can 
still be safely employed, since the probability distribu- 
tion P{E) is appreciably different from zero only in a 
support ~ [0 : i?c], for all values of g [see Fig. Q], and 
the two singularities in the densities of states and i^t, 
located at Ag and Ax respectively, overlap only through 
the high-energy tail of P{E), without significant contri- 



butions to the integral for Fqp. We then get 

d£:d£:' vs{E)r'"^{E)vT{E')P{E,E') 



As J A 



dE ,ys{E)r^{E)- / dE',,T{E')P{As,E') 



oc n 



neq 



/ dE' i^T {E')P {As, E') 

J At 



(38) 



where we used the notation P{E,E') = C{E,E')P{E - 
E'). The resulting rqp((7, Ag, At) is shown as a func- 
tion of g for some specific choices of Ag and Ax in 
Fig. [5] We choose to plot the dimensionless quantity 
Fqp = hTqp/{gTAs), meaning that the quasiparticle tun- 
nel rate is measured in units of As/h, and has been di- 
vided by gx- A superconducting gap of 2 K corresponds 
to a frequency of 40 GHz, and for gj- = 10^ the values 



shown in the figure indicate then Fq 



10 MHz. 
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FIG. 5: Dimensionless quasiparticle tunneling rate as a func- 
tion of environmental dimensionless conductance g. The red 
curve refers to the case As = 10-Ec, At = 5Ec, and no ob- 
servable dependence on g is noticed at this scale. The cyan 
curve has been obtained for As — 100-Ec, At ~ 50Ec, and 
the corresponding values of Fqp are slightly lower in this case. 

In the more interesting case Eq ^ SA (realized for ex- 
ample for Ag ~ 100 /xeV and C ~ IfF) the environment 
can couple energy regions where the singularities in the 
density of states of the two subsystems occur. Now the 
approximations adopted above are not justified anymore, 
especially for small values of g, and one needs in principle 
to solve the full two-dimensional integral in Eq. [25] The 
problem then is, without the decoupling of the integrals 
we cannot extract anymore a factor n"""^. We then need 
an explicit estimate for the unknown term /g'^'^. This can 
be done by assuming that the quasiparticles, while still 
being in thermal equilibrium at temperature T, are out 
of electro-chemical equilibrium, and the excess quasipar- 
ticle density can be accounted for by an effective 
chemical potential shift fl: 

- / di? ^g(g) r , .1 • (39) 
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To lowest order in temperature, we can connect fi directly 
to n'""i afP 



/i ~ Tin 



(40) 




FIG. 6: Same as in Fig. [S] but with different parameter 
values. The cyan curve corresponds to As = 2Ec, At = Ec 
and the red curve refers to the case As — Ec, At ~ 0.5Ec- In 
the first case, (As — At) equals Ec and leads to an unbounded 
increase in Fqp for gr — >■ 0. In the second case (and in general 
for Ec > 5 A) one observes rqp((7 — >■ 0) — >■ because the 
energy exchange Ec provided by the environment is too large 
to be absorbed by SA. 

Inserting the calculated fl in the formula for Fqp and 
performing the double integration, one can get numeri- 
cal estimations for any value of the parameters Ag/Ec 
and At/Ec- In Fig. [6] we report (cyan curve) the re- 
sults for the "worst" case (As — At) = Eq- One can see 
that in the limit 5 — >■ the quasiparticle poisoning rate 
is strongly enhanced, due to the perfect coupling of the 
two singularities in the density of states. However, this 
regime is difficult to attain, and the strong increase in 
Fqp is localized at g ~ which requires unrealistic envi- 
ronmental resistances R 3> Rq. In conclusion then, this 
issue should not represent a problem. 

In the regime (As — At) < Ec (red curve), the envi- 
ronment provides for g "^2 (see Fig. |4]) a typical energy 
larger than the "energy distance" between the two sub- 
systems, and since vt{E < At) = 0, smaller values for 
Fqp are obtained for decreasing g. In the limit g — > 
we have P{E) oc 6{E — Eq) and the result of integra- 
tion is suppressed to zero. Note that for 5 — > 00 the 
two curves of Fig. [6] approach each other, because in that 
limit P{E) is peaked \n E = and the exact position of 
At with respect to As becomes irrelevant. 



IV. QUANTITATIVE CONSIDERATIONS 

The final estimations strongly depend on the value of 
the tunneling resistance i?T which enters the expression 



for the poisoning rate. As anticipated above, such val- 
ues are different for different experiments, ranging from 
10 to ^ 10"' Vl. By looking at the expression for 
the quasiparticle tunnel rate, Eqs. (17) and (34) , one 



could conclude that large tunnel resistances (low g^) are 
desirable so that Fqp is reduced. But as we already com- 
mented, by the same token also Cooper-pair tunneling 
would be suppressed, and hence the proximity-induced 
gap would get reduced. Analytical calculation^ have 
shown that the pairing potential amplitude Apr induced 
in the proximized system, in terms of the parental pair- 
ing amplitude As and of the microscopic tunneling rate 
Fo, is given by 



Apr 



-A 



s • 



To + As 

The tunneling rate for bare electrons is evaluated as 

Fo = -^w^m , 



(41) 



(42) 



so that, using the definition Eq. (15) for i?T, one can 
relate Fq and i?T as: 



To 



(43) 



87ri?T i't(O) ' 

In the low transparency limit, Fq ^ As, the proximity 
gap is set by Fq, see Eq. (41), and is therefore rather 
small. On top of that, the topological gap is further 
reduced due to the Rashba and Zeeman interaction: 

At = ^==^==Apr = ^==Apr , (44) 



X 



with X = Vz/ioLk-p) quantifying the ratio between Zee- 
man splitting and typical spin-orbit interaction. Note 
that one always has At < Apr < As- Then, assuming 
the most favorable situation afcp ^ Vz (not so easy to 
achieve experimentally yelP^ and thus At — Apr, the 
requirement of a minimum topological gap of 100 mK 
translates into the condition Fg ~ Apr ~ 100 mK. 

As a final step, we need to estimate the normal-state 
density of states i^t(O) '^"^ the topological wire. To do so, 
we use the fact that the desired chemical potential has 
to lie in-between the gap opened by the Zeeman interac- 
tion added to the Rashba helical bands (at least in the 
simplest, ideal one-channel model). Using the dispersion 
relation 



e±{k) = rfc72m± 



U2 



(45) 



and requiring that the chemical potential lies in the mid- 
dle of the Zeeman gap, for instance halfway between 
£_(0) and £+(0) (as shown with red points in Fig. [7|, 
one gets the simple condition = 0. The ID density of 
states per unit volume at this energy is 



P^(e = 0) 



de_(fc)/d/c 
2 



e=0 



(46) 



k—ko 
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FIG. 7: Dispersion relation in a one-dimensional wire in 
the presence of Rashba spin-orbit and Zeeman interaction. 
The gap at fc = is entirely due to the Zeeman energy Vz. 
For afcp 2> Vz the position of the two minima e = eo is ap- 
proximately given by ±fcso = am/ti^ . The topological regime 
requires having the chemical potential lying inside the gap, 
as shown here; fco denotes the point at which the dispersion 
crosses the representative mid-gap level e = 0. 



where fco satisfies £_(A;o) = 0, see Fig. [t] Insertion of the 
expression for fcg in Eq. ( 46 1 leads to 



P^(e = 0) 



72 ^i + v/T+x^ 



(47) 



with X defined above. In the considered limit afcp 3> Vz 
and thus x <C 1 the density of states per unit volume is 
approximately given by 



i>^(e = 0) 



(48) 



The spin-orbit interaction strength a ranges from 0.00075 
eV-A in GaAs quantum well^ to 0.1 eV-A in In- 
GaAs quantum wells--?, or even more in heavier-element- 
wires such as InSbSl. We can therefore conclude that 
1^.5^(0) in the simple one-channel case varies between 10 
and 10^ (/im-K)~^. We choose the average value of 
~ 10^ (/im-K)~^ and a typical wire length of 1 /LtnPi. 
By imposing the constrain Fq ~ 100 mK derived above, 
we obtain via Eq. ( 43 1 the final estimate for the tunnel 
resistance 



100 . 



(49) 



As calculated in the former sections this value corre- 
sponds to a quasiparticle tunnel rate of Fqp ~ 1—10 MHz, 
i.e. poisoning times of the order of /xs or less, which has 
to be compared with the typical time required for adia- 
batic qubit manipulation. The natural time scale which 
identifies the adiabatic regime is provided by the inverse 
topological gap 



Tad 



h 

At 



(50) 



which for the aforementioned reasonable estimate of min- 
imum gap At = 100 mK takes the value Tad — 1 ns. The 
requirement of adiabatic computation is then satisfied 
if operations are performed on a time scale Tcomp much 
longer than Tad. In turn, quasiparticle poisoning events 
must be rare events during the time of computation: 



Tad <^ Tc, 



< Ta 



(51) 



where we have introduced for convenience the quasipar- 
ticle poisoning time Tqp = 1/Fqp. Assuming an order of 
magnitude difference between successive time scales, the 



above condition Eq. (51) sets the upper limit for Fqp to 



10 MHz, which is in the range of values we found in our 
calculations for an average situation. This shows again 
that the phenomenon of quasiparticle poisoning is not 
at all marginal, and its relevance should be assessed case 
by case. 

For example, for the only experimental results avail- 
able so far (Ref. [31]), the proximity effect is not very 
effective and the observed proximity gap is about one 
tenth of the bulk superconducting gap (which is however 
large in this case). On top of that, the spin-orbit energy 
is much smaller than the Zeeman energy in the topolog- 
ical phase, reducing the topological gap by an additional 
factor (approximately a factor 5 at the onset of the topo- 
logical transition). 

Note that Tad is set by the value of At, whereas Tqp is 
ultimately determined by Apr (via Fq and Rt) and does 
not depend on the physical properties of the topological 
nanowire (except for the density of states contained in 
-Rt)- Hence, the parameter regime akp ^ Vz is less 
favorable, not only due to the fact alone that one gets 
smaller values of the topological gap, but also because 
the adiabatic time scale is increased while the poisoning 
time remains constant. 

Working in the multi-channel regime would even be 
less favorable, since the density of states i?^ in the wire 
would be noticeably increased, and to maintain the same 
Fq the tunnel resistance R-j should be further decreased. 
A larger value of a would instead help in this direction, 
since it lowers (beyond increasing the topological 
gap). 

Also in the opposite limit of a transparent interface, 
Fq 3> As , where the proximity gap is essentially given by 
Ag, decreasing the quasiparticle tunnel rate is difficult. 
Equation ( 43 1 tells us again that for Fq > Ag ~ 1 meV, 



in order to suppress the factor Rq/Rt one would need 
unrealistically low values of the wire density. 

A possible improvement could be provided by the finite 
charging energy of the nanowire, which raises the energy 
of all the states and lifts the huge degeneracy of quasi- 
particle states close to At- For a single pair of Majorana 
states, the charging energy also introduces an undesired 
splitting between the filled and unfilled zero-energy state. 
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But one can then work with two wires and four Majorana 
states, two of which remain degenerate even in the pres- 
ence of a charging energjf^H. 

V. CONCLUSIONS 

In summary, we have calculated the tunnel rate ( "poi- 
soning") of quasiparticles from a bulk superconducting 
reservoir to a semiconducting nanowire, which becomes 
also superconducting due to proximity effect. Under ap- 
propriate conditions, the nanowire is in a topological su- 
perconducting state, hosting a Majorana state at each 
of its ends, which could be used for topological compu- 
tation. Using quantitative results from recent experi- 
ments on the density of excess quasiparticle in supercon- 
ductors, we have shown that the poisoning of the wire 



could represent a serious problem, with Majorana-qubit 
lifetimes which range from 10 ns to 0.1 ms, depending 
on many physical parameters. Since some of these pa- 
rameters cannot simply be adjusted independently, find- 
ing a suitable configuration which minimizes the poison- 
ing phenomenon requires a fine-tuning of the coupled 
nanowire-superconducting system more delicate than one 
could have expected. 



VI. ACKNOWLEDGMENTS 

We thank Luka Trifunovic for useful help with the nu- 
merical calculations. This work has been supported by 
the Swiss SNF, NCCR Nanoscience, NCCR QSIT, and 
the EU project SOLID. 



1 A.Y. Kitaev, Phys. Usp. 44, 131 (2001). 

^ M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 
103, 020401 (2009); M. Sato, Y. Takahashi, and S. Fuji- 
moto, Phys. Rev. B 82, 134521 (2010). 

^ R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010); Y. Greg, G. Refael, and F. von 
Gppen, Phys. Rev. Lett. 105, 177002 (2010). 

^ J. Alicea, Phys. Rev. B 81, 125318 (2010). 

^ L. Fu and G.L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

® A.R. Akhmerov, Phys. Rev. B 82, 020509(R) (2010). 

^ G. Goldstein and C. Ghamon, Phys. Rev. B 84, 205109 
(2011). 

® J.G. Budich, S. Walter, and B. Trauzettel, Phys. Rev. B 

85, 121405(R) (2012). 
^ B. van Heck, F. Hassler, A.R. Akhmerov, and C.W.J. 

Beenakker, Phys. Rev. B 84, 180502(R) (2011). 
^° L. Fu and G.L. Kane, Phys. Rev. B 79, 161408(R) (2009). 

F. Hassler, A.R. Akhmerov and G.W.J. Beenakker, New 

J. Phys. 13, 095004 (2011). 

P.J. de Visser, J. J. A. Baselmans, P. Diener, S.J.G. Yates, 
A. Endo, and T.M. Klapwijk, Phys. Rev. Lett. 106, 167004 
(2011). 

J. Aumentado, M.W. Keller, J.M. Martinis, and M.H. De- 
voret, Phys. Rev. Lett. 92, 066802 (2004). 
A.J. Ferguson, N.A. Gourt, F.E. Hudson, and R.G. Glark, 
Phys. Rev. Lett. 97, 106603 (2006). 

M.D. Shaw, R.M. Lutchyn, P. Delsing, and P.M. Echter- 
nach, Phys. Rev. B 78, 024503 (2008). 
M. Zgirski, L. Bretheau, Q. Le Masne, H. Pothier, D. Es- 
teve, and G. Urbina, Phys. Rev. Lett. 106, 257003 (2011). 

^"^ L. Sun, L. DiGarlo, M.D. Reed, G. Gatelani, L.S. 
Bishop, D.I. Schuster, B.R. Johnson, G.A. Yang, L. Frun- 
zio, L.I. Glazman, M.H. Devoret, and R.J. Schoelkopf, 
|arXiv :1112.2621 

^® M. Tinkham, Introduction to Superconductivity (McGraw- 
Hill, New York, 1996), 2nd ed. 

1^ S.B. Kaplan, G.G. Ghi, D.N. Langenberg, J.J. Ghang, S. 
Jafarey, and D.J. Scalapino, Phys. Rev. B 14, 4854 (1976). 
J.D. Sau, R.M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. B 82, 094522 (2010); J.D. Sau, S. Tewari, R.M. 
Lutchyn, T.D. Stanescu, and S. Das Sarma, Phys. Rev. 



B 82, 214509. (2010); J.D. Sau, S. Tewari, and S. Das 
Sarma, Phys. Rev. B 85, 064512 (2012); T.D. Stanescu, 
R.M. Lutchyn, and S. Das Sarma, Phys. Rev. B 84, 144522 
(2011). 

S. Gangadharaiah, B. Braunecker, P. Simon, and D. Loss, 
Phys. Rev. Lett. 107, 036801(2011); E.M. Stoudenmire, 
J. Alicea, G.A. Starykh, and M.P.A. Fisher, Phys. Rev. B 
84, 014503 (2011). 

J. Ahcea, Y. Greg, G. Refael, F. von Gppen, and M.P.A. 
Fisher, Nat. Phys. 7, 412 (2011). 

G.-L. Ingold and Yu.V. Nazarov, Single Charge Tunneling, 

edited by H. Grabert and M.H. Devoret, NATO ASI Series 

B Vol. 294 (Plenum, New York, 1992), p. 21. 

R. M. Lutchyn, Ph.D. thesis, Yale Univ. 

R. Lutchyn, L. Glazman, and A. Larkin, Phys. Rev. B 72, 

014517 (2005). 

'^^ G. Gatelani, J. Koch, L. Frunzio, R.J. Schoelkopf, M.H. 
Devoret, and L.I. Glazman, Phys. Rev. Lett. 106, 077002 
(2011). 

G. Gatelani, R.J. Schoelkopf, M.H. Devoret, and L.I. Glaz- 
man, Phys. Rev. B 84, 064517 (2011). 
J.M. Martinis, M. Ansmann, and J. Aumentado, Phys. 
Rev. Lett. 103, 097002 (2009). 

2^ M. Hofheinz, E.M. Weig, M. Ansmann, R.G. Bialczak, E. 
Lucero, M. Neeley, A.D. O'Gonnell, H. Wang, J.M. Marti- 
nis, and A.N. Cleland, Nature (London) 454, 310 (2008). 

^° O. Naaman and J. Aumentado, Phys. Rev. B 73, 172504 
(2006). 

J.T. Pehonen, J.T. Muhonen, M. Meschke, N.B. Kopnin, 
and J.P. Pekola, Phys. Rev. B 84, 220502(R) (2011). 
A.D. Gorcoles, J.M. Ghow, J.M. Gambetta, G. Rigetti, 
J.R. Rozen, G.A. Keefe, M.B. RothweU, M.B. Ketchen, 
and M. Steffen, Appl. Phys. Lett. 99, 181906 (2011). 
J.M. Martinis, M. Ansmann, and J. Aumentado, EPAPS 
Document No. E-PRLTAO-103-051936. 
V. Mourik, K. Zuo, S.M. Frolov, S.R. Plissard, E.P.A.M. 
Bakkers, and L.P. Kouwenhoven, Science 336, 1003 (2012). 
L. Meier, G. Sahs, I. Shorubalko, E. Gini, S. Schon, and 
K. Ensslin, Nat. Phys. 3, 650 (2007). 

J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 
Rev. Lett. 78, 1335 (1997). 



11 



J.D. Sau, B.I. Halperin, K. Flcnsberg, and S. Das Sarma, 
Phys. Rev. B 84, 144509 (2011). 

For p— wave pairing the u and v amplitudes become 2x2 



matrices, but in our effective spinless case we assume we 
can still use scalar s— wave-like coefficients. 



